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Abstract. In this second paper, in a series of two, we determine the properties of the stellar torus that we showed 
in the first paper to result as a product of two merging black holes. If the surrounding stellar cluster is as massive 
as the binary black hole, the torque acting on the stars ejects a fraction which extracts the binary's angular 
momentum. After the black holes coalesced on scales of ~ 10 7 yr, a geometrically thick torus remained. In the 
present article we show that a certain fraction of the stars has winds, shaped into elongated tails by the central 
radiation pressure, which are optically thick for line of sights aligned with them. These stars are sufficiently 
numerous to achieve a covering factor of 1, so that the complete torus is optically thick. This patchy structured 
torus is then compared with observations. We find the parameters of such a torus to be just in the right range in 
order to explain the observed large column densities in AGN and their temporal variations on time scales of about 
a decade. Within this model the broad absorption line quasars can be interpreted as quasars seen at intermediate 
inclination angles, with the line of sight grazing the edge of the torus. The half-opening angle of the torus is wider 
for major mergers and thus correlates with the central luminosity, as has been suggested previously. In this picture 
the spin of the merged black hole is possibly dominated by the orbital angular momentum of the binary. Thus 
the spin of the merged black hole points into a new direction, and consequently the jet experiences a spin-flip 
according to the spin-paradigm. This re-orientation could be an explanation for the X-shaped radio galaxies, and 
the advancing of a new jet through the ambient medium for Compact Symmetric Objects. 

Key words. Black hole physics - Stellar dynamics - Galaxies: active - Galaxies: interactions - Galaxies: jets - 
Galaxies: nuclei 



■ 1. Introduction gions of the accretion disk, which produces the optical- UV 

emission in the outer parts, extending to about 10 pc. 

Active galactic nuclei (AGN) produce very high luminosi- Qn scales up to g.Olpc the broad emission line region 

ties in a very restricted volume, probably via an accretion (BLR) is found (Peterson 1993) with line widths corre- 

disk around a massive black hole (BH) in which shear sponding to doppler-broadening with velocities of 1500 up 

stresses cause matter to sink down in the potential of to 30 000 km/s, probably caused by line emission in clouds 

the BH. These luminous nuclei appear in many different moving at these speeds while these lines are varying on 

flavours what the unification scheme traces back to a non- scaleg of weeks Qr montnSj in agreement with the exten- 

spherical symmetry of the nucleus, which is spatially not gion of the BLR the narrQW emission lines do not show 

resolved. Thus the same object looks different if viewed such variat ions. Their line widths imply velocities of the 

from different directions, leading to various classifications. order of 10Q0 km/s and hence place the narrow emission 

Such an AGN we tried to illustrate in Fig. 1 in une region (NLR) at distances up to the kilo parsec scale, 
double logarithmic scales. In the very center is a mas- 
sive BH of about 10 s M with a Schwarzschild radius of 

about 10~ 5 pc. The observed rapid variation of X-rays In the infrared (IR) the spectrum shows a bump be- 
on scales of 10 4 s restricts their origin to a volume of tween 2 jum and 1 mm which is thought to be due to re- 
about 3 x 10 12 m radial extension, which corresponds to radiation by dust in a distance of some parsec from the 
the last stable orbit of a black hole of — 3 x 10 8 M Q . center (Chini et al. 1989; Sanders et al. 1989) and which 
The X-rays are thought to be generated near the inner re- is distributed in a toroidal volume around the nucleus, 
whose symmetry axis is aligned with that of the accretion 
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Fig. 1. A sketch of the cylindrically symmetric AGN according to the description in the introduction and to the picture of a 
patchy torus in AGN, developed in paper I and this work, is shown. The cut shows the r — a-plane, both axes logarithmically 
scaled to 1 pc. The basic constituents are the central BH with a surrounding accretion disk, the jet perpendicular to the disk 
and the torus encircling this configuration. The dark patches in the torus indicate the clouds, made up by stellar winds in our 
proposed concept. The locations of the BLR, ESR and NLR are shown on the left, as well as the appearance of the nucleus as 
a function of the angle of the LOS to the jet-axis (on the right). See text for details. Note that a strictly spherical distribution 
is almost a square box in this diagram. 



The observed AGN can basically be divided in two 
different types, with Type 1 objects exhibiting all the fea- 
tures described above. Therefore in these objects the cen- 
tral parts are directly seen without absorbing matter in 
the line of sight (LOS). On the other hand in sources of 
Type 2, there is no convincing evidence for the directly 
visible optical-UV emission and the soft X-rays seem to 
be heavily absorbed (Lawrence & Elvis 1982; Mushotzky 
1982). No broad emission lines are detected. The unifica- 
tion scheme relates the different types to different orien- 
tations rather than intrinsic differences: The Type 2 clas- 
sified nuclei are seen at large inclination angles and conse- 
quently the LOS is blocked by the torus so that only the 
NLR is freely visible and in fact indistinguishable from 
that in Type 1 objects. 

Evidence for such obscuring tori is deduced from ob- 
served sharply defined conical/biconical nebulae which are 
suggestive of shadowing (Antonucci (1993) and references 
therein). The deficit of observed ionizing photons relative 
to what is needed to ionize the nebulae (Kinney et al. 
1991) evidently supports the idea of matter obscuring the 
central source, which is not obscured at the aspect angles 



where the ionized nebulae are located. Hence the ionizing 
photons escape from the center in a cone so that the op- 
tical, UV and probably soft X-ray radiation is collimated 
by the torus, escaping anisotropically through the polar 
cap regions. 

Even stronger support comes from spectrapolarime- 
try, which reveals in the polarized spectrum of Type 2 
sources a clear spectrum of Type 1 (Antonucci & Miller 
1985; Miller ct al. 1991). The power-law continuum has the 
same degree of polarization as the permitted broad lines, 
which also have normal equivalent width. The degree of 
polarization is found to be almost wavelength indepen- 
dent. This is interpreted as the reflection of the obscured 
Type 1 continuum by free electrons in the opening of the 
torus (Antonucci (1993) and references therein). This elec- 
tron scattering region (ESR) is located somewhere near to 
the BLR and extends to scales of the height of the torus 
(Taniguchi & Anabuki 1999), see Fig. 1. 

Thus for a small viewing angle with the source be- 
ing seen almost face on all features of a Type 1 spectrum 
are visible. As the viewing angle increases the LOS finally 
grazes the edge of the torus and the spectrum becomes 
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reddened. The central regions gradually become more ob- 
scured by the torus till the broad lines disappear and the 
source is classified as Type 2 when it is seen near edge on. 

As indicated in Fig. 1 , there are more criteria than the 
orientation to distinguish between the different objects, 
such as the luminosity or the radio properties, with radio 
loud galaxies having a strong jet streaming along the sym- 
metry axis. For a much more detailed explanation of the 
unified scheme see for example the reviews by Antonucci 
(1993), Urry & Padovani (1995) and Wills (1999). 

The torus is a very important component in this 
scheme. To support its vertical thickness previous models 
either use magnetic pressure (Lovelace et al. 1998) or ra- 
diation pressure (Pier & Krolik 1992b), which faces some 
problems for low luminosity central sources and smoothly 
distributed dust. Here we want to further develop our idea, 
introduced in a previous paper (Zier & Biermann 2001), 
that this torus is comprised of the winds of stars moving 
in the potential of a merging binary black hole (BBH) , al- 
lowing them to maintain the geometrical thickness of the 
torus. 

In that article we simulated the evolution of a stellar 
cluster under the influence of a massive binary black hole 
in its center, based on the two widely accepted assump- 
tions that (a) galaxies with an active nucleus harbour a 
supermassive black hole in their center and (b) that galax- 
ies frequently merge. The results clearly showed that the 
binary merges on scales of some 10 7 yr due to ejection of a 
fraction of stars, if the surrounding stellar cluster is about 
as massive as the binary, and ensuing emission of gravita- 
tional radiation. The stars which remained bound indeed 
are found to be distributed geometrically in the volume 
of a thick torus. In this second paper we will investigate, 
wether such a stellar torus can explain the presence of the 
ubiquitous torus surrounding apparently all AGN. 

The next section will briefly review the constraints put 
on the properties of the torus by observations. 

2. Properties of the torus 

Although the current telescopes are not able to spatially 
resolve the torus in AGN, conclusions could be drawn from 
observations in various wave bands about the torus' geom- 
etry as well as the matter it is made of. 

The inner radius 

Absorption of the primary optical-UV radiation by dust, 
which reemits it in the IR, heats the dust to its evapora- 
tion temperature (<~ 1500 K). Inside the evaporation ra- 
dius dust can not exist and Lawrence (1991) assumed it to 
be the inner radius of the dusty torus. This radius is about 
lpc, corresponding to the distance where the BBH be- 
comes hard (Milosavljevic & Merritt 2001) and what also 
marks the inner radius of the torus that we obtained in our 
simulations (paper I). This value is also in agreement with 
the inner radius obtained by Krolik & Begclman (1988) 



from the balance between cloud evaporation by central ra- 
diation and inflow by dissipative processes in the torus. In 
several Type 1 AGN variations in the near-infrared (NIR) 
emission have been observed to lag those at shorter wave- 
lengths (Clavel et al. 1989; Sitko 1991; Baribaud et al. 

1992) . These time delays place the IR emission source in 
a distance from the central engine which is consistent with 
that of optically thin nuclear heated dust at its evapora- 
tion temperature. 

To explain the bump in the near-infrared in terms of 
thermal dust emission, a range in the grain temperature 
is required (1500 K < Td us t ^ 30 K), since isothermal 
dust close to its evaporation temperature yields a narrower 
bump than that observed (Sanders et al. 1989; Haas et al. 
2000). Thus the torus' temperature has to decrease from 
the inner edge to its outer regions (Niemeyer & Biermann 

1993) , but a very thick and compact torus (Pier & Krolik 
1992a, 1993) is not able to provide such a broad tempera- 
ture range. Instead a clumpy torus, as suggested by Krolik 
& Begelman (1988), has the advantage that dust in clouds 
can survive more easily the strong radiation field, and such 
a structure of the torus tends to increase the dust temper- 
ature in the outer parts (Efstathiou & Rowan-Robinson 
1995) where the clouds are exposed to central radiation 
through gaps in the inner parts of the patchy torus. This 
seems also most promising to Haas et al. (2000) in order 
to explain the spectral energy distribution (SED). 

Column density Nn 

On the basis of X-ray observation lots of evidence has 
been accumulated that the column density is in the range 
10 22 - 24 cm- 2 (Mulchaey et al. (1992) for UV-detected 
Seyfert 2s), or l-8x 10 24 cm" 2 (Krolik & Begelman 1988). 
A more recent study of the X-ray observations of a sample 
of Seyfert 2 galaxies in the 0.1 — 100 keV spectral range 
with BeppoSAX has been performed by Maiolino et al. 
(1998). The sources they studied were selected according 
to their [OIII] optical emission line flux after correction 
for the extinction deduced from the Balmer decrements. 
Hence this sample is thought to be fairly isotropic and 
avoids a bias against heavily obscured nuclei, as is proba- 
bly the case in former samples. Maiolino et al. (1998) find 
all the sources to be absorbed by column densities larger 
than 4 x 10 23 cm~ 2 and most even appear to be thick to 
Compton scattering with A^h > 10 25 cm~ 2 . Hence their 
results provide further and strong evidence for the unifi- 
cation scheme and indicate that the obscuration in Type 2 
AGN is probably much higher than deduced from earlier 
X-ray surveys. 

Such high column densities are also deduced from radio 
data. Conway & Blanco (1995) ascribe the absorption of 
the broad HI 21 cm line towards the compact radio nucleus 
of the FR II radio galaxy Cyg A to the obscuration by a 
torus in about 10 pc distance to a central BH of <~ 10 s M Q . 
Assuming a spin temperature in the range 8000—16 000 K 
the authors obtain a column density in the range 2 — 4 x 
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10 23 cm 2 , in very good agreement with the column 3.75 x 
10 23 cm -2 deduced from previous X-ray spectroscopy by 
Ueno et al. (1994). 

Both, X-ray and radio data show also evidence for a 
patchy structured torus (Boiler et al. 2002; Krichbaum 
et al. 1998), in agreement with the explanations for the 
IR range of the previous section. 

Opening angle 

The observed sharply defined conical/biconical nebulae in 
some AGN are suggestive of shadowing (Antonucci (1993) 
and references therein). These ionization-cones, being lo- 
cated in the opening of the obscuring torus and aligned 
with the inner radio jets (Nagar & Wilson 1999), repre- 
sent a lower limit of the half opening angle of torus in the 
range 35° to 50° (Pogge 1989; Wilson et al. 1993; Wilson 
& Tsvetanov 1994; Macchetto et al. 1994). 

Using the statistics of the observed number ratios of 
Seyfert 1 and 2s (Osterbrock & Shaw 1988; Huchra & Burg 
1992; Barthcl 1989; Willott et al. 1999) a half-opening 
angle of the torus in the range 30° — 60° is obtained. This 
matches very well the angles of the ionization-cones inside 
the opening of the torus. 

Covering factor 

According to the opening angle of the torus, it screens 
a considerable fraction of the solid angle, seen from the 
central source, and is geometrically thick. Blocking the 
view to the nucleus in Type 2 AGN, the torus, if comprised 
of individual clouds, must be optically thick for most of the 
lines of sight penetrating it. Thus the clouds constituting 
the torus must achieve a covering factor of about one, if a 
single cloud is optically thick, and more otherwise. 

Hence the observations demand a clumpy-structured 
torus with an inner radius of about 1 pc, that is optically 
and geometrically thick with a half-opening angle of ~ 45° 
and therefore has a covering factor of order unity. 



3. Kinematics and density distribution of the stars 

Our simulation of a stellar distribution in the potential of 
the binary (discussed in detail in paper I, Zier & Biermann 
(2001)) has been performed for three different mass-ratios 
of the BHs (q = 1, 10, 100), with 

Mi 

9 = w a > 

where Mi > M 2 , always. While we focused on the ejected 
fraction of the stars in that article we will discuss in this 
section the bound population. For clarity we first briefly 
summarize just a few results from paper I, using a mass- 
ratio of the BHs of q = 1. 
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Fig. 2. The stellar number density distribution is shown at the 
beginning (t = 0) and end (t = t m ax) of the simulation. Finally, 
at fmax, the bound star distribution follows a powerlaw with 
index —4 in the heated region (f > 50), while an index ~ —2 
is maintained in the range 10 < f < 50. The inner parts are 
scoured out due to ejection of stars by the binary, and a max- 
imum emerges at f ~ 3, showing the torus-like configuration 
the stars assume. For r < 0.7 a cusp of stars bound to Mi 
only is left, f is normalized to 1 pc and h such that the area 
underneath the solid line is 1. 

3.1. Dynamics of the bound stars 

In order to write the equations in their dimensionless form, 
denoted by a tilde " ~ " on top of the quantities, we used 
the following normalization parameters: 

ro = a = 1 pc , 

, _ / r- 3 
\/G(Mi+M 2 )' 

ma 2 



The semi-major axis of the BBH is denoted by o, and for 
the mass of a typical star m we used one solar mass. 

During the simulation the initial singular isothermal 
sphere evolves finally into a distribution which is depleted 
in the inner regions (see Fig. 2). The stellar density peak 
at about f = 0.5 is made up by a fraction less than 0.4% 
of the final distribution and is bound to the primary black 
hole only. Since we expect this cusp to be much shallower 
if we would have taken into account the shrinking of the 
binary in the calculations, we will not address it further 
in our discussion. 

After a gap in a distance corresponding to the radial 
distance of the secondary black hole to the center, the 
density increases again and peaks at f <~ 3. In the range 
10 ^ ^ ^ 50 the density follows the initial powerlaw of 
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Fig. 3. The initial conditions (t = 0) of the star cluster have 

been set so that the velocity anisotropy /3 va = 1 -^tftt^- 

is neither radially (> 1) nor tangentially (< 1) anisotropic. At 
t — tmax the velocity anisotropy of the bound stars is a strongly 
increasing function of the radius. In the range r < 50 they are 
tangentially anisotropic, while in the heated region (r > 50) 
the stars are moving on radially anisotropic orbits. 



the isothermal sphere (h oc f~ 2 ) and for radii bigger than 
the outer limit of the initial distribution (f > 50) approx- 
imates a powerlaw with index —4. 

In paper I we showed that the stars with small angular 
momenta, i.e. small pericenters, get close to the orbit of 
the secondary black hole where they are strongly inter- 
acting with both BHs, so that they become ejected. Only 
stars moving on orbits with sufficiently large pericenters 
(f_ > 2) can remain in the central regions. The smaller 
their radii, the more they are tangentially anisotropic (see 
Fig. 3) and the less eccentric are their orbits. A fraction of 
the stars which are interacting with the BHs will be heated 
to larger distances and stays bound to the binary instead 
of becoming ejected. Such stars did not gain as much en- 
ergy as the ejected ones and contribute to the distribution 
at distances larger than the extension of the initial distri- 
bution, which terminates at f = 50 (see Fig. 2). 

In addition to these stars we will find in this region 
(f > 50) also those which happen to be close to their 
pericenters at the very beginning of the simulation due 
to the choice of the initial conditions. As time proceeds 
they are moving on their orbits without being noticeably 
affected by the secondary BH, whose perturbation to a 
point-mass potential is negligible at large distances. Thus, 
keeping their energy and angular momentum nearly con- 
served, these stars diffuse into the outer regions (f > 50). 
The expansion of the initial distribution is just a process 
of relaxation and would have occurred also in a point-mass 



Fig. 4. Plotting the angular momentum L versus the distance 
f at i max shows all stars to respect the upper limit indicated 
by the dotted line. Only some stars in the cusp region (r < 1) 
have larger angular momenta, because they are bound to one 
BH only, with the other perturbing their orbits, so that the 
assumption of a point-mass potential does not apply at these 
small distances. For f > 50 two populations can be distin- 
guished, one clustering around L m 2 (heated stars) and the 
other having an upper limit of L « 10, fading towards smaller 
values (relaxation process). See text for details. 



potential. Thus the region at r > 50 is partly populated 
by stars which did not interact with the binary. 

Both these populations can be clearly distinguished in 
Fig. 4, where we plotted the angular momentum versus the 
radius for each single star at the end of the simulation. 
For f > 50 one population is clustering around L ~ 2 
and the other having an upper limit at L w 10, fading 
towards smaller angular momenta. Both are extending to 
radii of about 10 4 , with the former population being more 
pronounced at such big distances. 

In a point-mass potential the energy reads in dimen- 
sionless units 



1 



(1) 



with the maximum energy of a bound star just below zero. 
To maximize the angular momentum, all its kinetic energy 
is assumed to be stored in tangential motion so that v r = 
0. Solving now for the angular momentum gives 



L < V2r 



(2) 



as the upper limit for bound stars. This power law with 
index —1/2 is indicated in Fig. 4 by the dashed line. 

However, the nature of the two populations at r > 50 
can be understood easily if we calculate the maximum 
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Fig. 5. Cuts through the stellar density in the comoving frame are shown with contours scaled logarithmically. The right panel 
displays the equatorial plane (BHs marked by black spots). Perpendicular to it the x — 0-plane is shown (left panel), with the 
j/-axis drawn as dashed line, so that the BHs are in front and behind the paper-plane. The initial distribution is a gaussian and 
the mass-ratio is 1. After initially stars close to the orbits in 0.5 pc distance from the center are ejected and the polar regions 
are depleted, finally a torus emerges in the equatorial plane at r ~ 3pc. An expanded version of this diagram for the mass-ratio 
10 is shown in paper I. 



possible angular momentum for the stars as function of 
their initial pcricenter. In dimcnsionless form this is 



(3) 



to time from the beginning till the end of the simulation 
(^max = 5x 10 5 ) this distance is 



1 



In order to obtain an upper limit for the angular mo- 
mentum of the diffused population, we assume that some 
stars at f = 50 arc initially at their pericenter and move 
on highly eccentric orbits (e — > 1). Solving then Eq. (3) 
for L yields an upper limit of 



(4) 



If we use for the initial pericenter r_ = 2 instead of 50 
we obtain Ly lm = 2. This clearly shows that the stars at 
r > 50, whose angular momenta tend to 10, initially have 
been close to f — 50 and thus diffused into this region. On 
the other hand the stars with angular momenta L sa 2 had 
initial pericenters close to the radius of the orbit of the sec- 
ondary BH. Here they strongly interacted with the BBH 
before they were heated to distances up to 10 4 (Fig. 4). 
The heated stars in fact are very similar to the ejected 
ones which have been discussed in detail in paper I. Their 
velocity distribution is also radially anisotropic (/3 va > 0, 
see Fig. 3) and they are moving on highly eccentric orbits. 
The heated stars just did not gain sufficient energy as to 
escape from the binary. 

This can be confirmed by computing the maximum 
distance a star can travel, when it is moving at escape 
velocity, t> osc = v/2/f. If we integrate v csc with respect 



3 - 



-3/2 

r. 



10 4 



(5) 



corresponding to the largest radii in Fig. 4. For the initial 
radius n we have inserted values ranging from to 50, but 
since the first term in the brackets dominates the expres- 
sion it does not make a difference whether the stars started 
in the center or from the edge of the initial distribution. 

The bound stars in the intermediate range eventually 
assume a torus-like structure as is shown in Fig. 5, which 
is taken from paper I. This plot displays cuts through the 
3-dimensional stellar density distribution in the comoving 
frame of binary after about 10 7 yr when the BHs have 
merged. The right panel shows the equatorial plane (z — 
0) with the BHs marked by the black points, M\ to the 
right. The slice in the left panel is perpendicular to the 
equatorial plane and contains the x = 0-plane, with the y- 
axis indicated by the dashed line. While the stars initially 
have been distributed according to a gaussian, the central 
parts have been scoured out and a torus in the equatorial 
plane of the binary with a radius of about 3pc is left. 
The basic topology of the final density distribution can be 
understood in physical terms, when we consider the torque 
exerted by the two black holes on an orbiting star. Fig. 3 
in paper I shows of the normalized torque the component 

L z ,i/r — r ifi sin 2 (9) relative to the z-axis of the BBH as a 
function of time for different angles 9 between the rotation 
axis of the BBH and the symmetry axis of the star's orbit. 
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Fig. 6. The mean of the squares of the velocity components in 
spherical coordinates are displayed at i max . In the range of the 
initial distribution (f < 50) both tangential components exceed 
the radial one, while at larger distances it is just the opposite, 
in accordance with Fig. 3. In the range from the inner edge of 
the torus at f — 1 to the maximum of the number density at 
f = 3 (Fig. 2), the azimuthal component is larger than (%}, 
in favour of the toroidal distribution of the stars. At larger 
distances these components become indistinguishable. 

The bigger the angle 9 is (i.e. for orbits through the polar 
regions), the more is the star's trajectory disturbed by the 
influence of the two BHs. The cumulative effect of these 
large excursions in the polar cap regions deplete the stellar 
population, leaving a torus behind. 

The toroidal structure which the bound population 
assumes once the BHs are merged is confirmed also by 
the kinematics of the stellar distribution. The tangential 
velocity components of a single star moving on an orbit 
which is inclined by the angle 9 to the equatorial plane 
are 

U0 = v t cos 9 V 
v e = v t sin V , 

where 9 V is the angle enclosed by the binary's rotation axis 
and the velocity-vector v of the star. Thus the total tan- 
gential velocity component is v t = (v^ + v^) 1 ^ 2 . For stars 
constrained to a torus-like volume, as shown in Fig. 5, we 
expect their azimuthal velocity component to exceed 
the polar component vg. This is verified in Fig. 6, where 
we have plotted the mean square of each velocity com- 
ponent versus the radius: At the inner edge of the torus 
between 1 and 2pc distance (Figs. 2 and 5) the mean 
square of the azimuthal component is a factor of about 10 
larger than (v" 2 }. With increasing radius both components 
approach each other and at about 8pc they are almost 
indisinguishable. 



4. Stellar winds comprising a patchy torus 

Now that eventually the bound stars indeed constitute a 
torus-like distribution which is geometrically thick, just as 
required by the unification scheme, we may ask about the 
effects of this structure on the radiation emitted from the 
center. Of course the stars alone are not able to obscure 
this radiation, since at a distance of about r — 3 pc a 
number of order of at least (r/i?*) 2 w 10 12 stars would 
be required to achieve a covering factor of order unity, 
using giant stars with a radius of 10 11 m. But if there is a 
sufficiently large number of stars in the torus with dense 
winds which absorb a significant fraction of the central 
radiation, the complete torus might be opaque. To get an 
estimate of the optical depth and density of the torus we 
have to examine its constituents, the stars and their winds, 
more closely. 

Solar type stars are not appropriate for the obsuration 
of the AGN, since their winds are too weak. But there 
are various other stars which can account for the obscu- 
ration and which we will call 'obscuring stars' (OS) in the 
following. Using evolutionary tracks, Young et al. (1977) 
showed that the total fraction on the giant branch in a 
cluster ranges between 1% and 2% for an interval 1 to 
5 billion years after coeval star formation. In a first pa- 
per Shull (1983) assumed a dense stellar cluster to consist 
of two components, one evolved component of red giants 
(1-M©, 100 R & ) comprising 1% of the stars, and a main- 
sequence component. In a later paper (Voit & Shull 1988) 
the fraction of supergiants in a thermal cluster was as- 
sumed to comprise only 0.01%, giants (1M , 10 i?©) 1% 
and the rest of the cluster to be made up by main-sequence 
stars. This would be only a small fraction with stars hav- 
ing mass- loss rates of typically 10 -5 M©/yr. On the other 
hand the observed FWHM within the inner 160 pc of both 
Ha and [N II] and the equivalent width of the Ca II lines 
are indicating the presence of an important population 
of red supergiants in the nucleus of NGC 6951 (Perez 
et al. 2000). The current understanding of stellar winds 
is based on observations in the solar neighbourhood and 
so far there are no results concerning the possible effect of 
the intense radiation of the AGN on the wind structure. 
Thus stars in AGN might be very different from stars in 
the outer galaxy (Alexander & Netzer 1994, 1997). There 
are theories that the conditions in AGN are likely to en- 
hance the mass-loss rate of irradiated stars near the cen- 
ter (Edwards 1980; Shull 1983; Voit & Shull 1988; Tout 
et al. 1989). While this radiation might not significantly 
increase the mass-loss rate (Voit & Shull 1988), Biermann 
(1989) and Stecker et al. (1991) suggested a significant 
neutrino flux in AGN which may transform low-mass stars 
into red giants (MacDonald et al. 1991). Since we consider 
the the torus to be a consequence of two merging galaxies 
with central massive black holes, besides stars we also ex- 
pect big amounts of gas to be brought to the center, and 
therefore an enhanced formation rate of young stars with 
high mass-loss rates. This is supported by Perez et al. 
(2000), who find evidence that star formation occurs in 
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Fig. 7. This sketch illustrates how the radiation pressure of a 
strong source, located in a large distance d (d ^> 7? w ,ex) left 
from the star, shapes the wind into a comet-like tail, pointing 
radially away from the source. The radius of the cross-section 
area of the wind is _R W , and i? WiOX denotes the radius of a freely 
expanding wind. 



bursts and continuously inwards to the nucleus from a 
5arcsec radius in NGC 6951, and Oliva et al. (1995), who 
investigated red supergiants as starburst tracers in galac- 
tic nuclei. They find a large value of Lr/M in NGC 1068, 
confirming the presence of a decaying starburst in the cen- 
tral regions (R < 3"). We will condense the previous ar- 
guments into the simple assumption, that the obscuring 
stars with strong winds amount to a fraction of 1% of the 
central cluster. 

The radiation pressure of a bright central source lim- 
its the wind-radius in its extension and will shape it into 
an elongated tail, pointing radially away from the cen- 
tral source. This is illustrated for one star in Fig. 7. The 
large light-grey shaded circular area represents the freely 
expanding wind of the star in its center. The dark-grey 
area shows the smaller extension of the wind, limited by 
a central source far left from the star, which shapes the 
wind into a comet-like tail. 

Such bending of the wind into a tail by radiation pres- 
sure can be deduced from a crude and simple estmate. 
According to Mathis et al. (1977), Graphite with its evap- 
oration temperature T cvap ~ 1500 K is the major contrib- 
utor to the extinction of the radiation outside the evapo- 
ration radius r cvap . This corresponds to the distance from 
the center, where the absorption of the central UV radia- 
tion by dust grains equals the rate at which it is reradiated 
as thermal IR radiation, such that the equilibrium temper- 
ature is the evaporation temperature. The absorption effi- 
ciency of the grains in the IR is much smaller than in the 
optical/UV, so that the reemission in the near IR (NIR) 
must be optically thin. Barvainis (1987) approximated the 
IR-absorption efficiency by a powerlaw {Q v cx v 1 ) with 
the index in the range 1 to 2. This yields a spectrum with 
a considerably narrower bump at 2 /im than is observed 
(Sanders et al. 1989; Haas et al. 2000) and which corre- 



sponds to the optically thin emission peak of the hottest 
grains at a temperature T gr ~ 1500 K. Hence, any model 
which seeks to explain the NIR bump in terms of thermal 
dust emission requires a range in the grain temperature. 
This is naturally included in our model, since the heated 
dust will be in equilibrium at different temperatures at 
different distances from the central radiation source. For 
the grain temperature as function of the distance r to the 
center Barvainis (1987) obtained 



T, r = 1650 



L 



uv,46 



e -5r uv /28 K> 



(6) 



with r uv as the optical depth of dust in the UV and L 



uv,46 



being the central UV luminosity in units of 10 erg/s. 
is the distance normalized to 1 pc. Using T evap for the 



' pc 



grain temperature and solving this equation for r gives 
the evaporation radius 



r - 4 05 T 1/2 T" 14/5 r>c 

'evap — ^ uv ,46 J evap,3 P L ' 



(7) 



where T emPj 3 denotes the evaporation temperature in 
units of 1000 K. 

Assuming a star to be located inside r evap where the 
Thomson cross-section cr th applies, its wind extension can 
be calculated using the same Ansatz as for the Eddington 
limit, i.e. equating the force of the central radiation acting 
on the wind particles with the kinetic force of the wind, 



^uvCth 

47rr 2 c 



(8) 



With typical values of red giants for the wind velocity i> w 
(Knapp & Morris 1985; Winters et al. 1994; Hoefner et al. 
1996) this equation can be solved for the extension of the 
wind, 



1.2 x 10" 



r 
2pc 



10 



km 



Vio 46 



pc 



(9) 



This is only a very rough estimate. In the following we 
are seeking for a more robust value of the wind radius, 
which is based on confirmed limits set by observations. As 
discussed in Sect. 2, these tell us that the torus is optically 
thick and the covering factor of its constituents, the stars, 
is about one. Therefore the number density of obscuring 
stars, integrated along the radial extension of the torus, 
yields a surface density of the number of stars (S os ) which 
has to equal the inverse of the cross-section area of the 
stellar winds, i.e. 



1 



(10) 



i? w denotes the radius of the cross-section area of the stel- 
lar wind perpendicular to its radial extension. To calculate 
the surface density we employ the number density pro- 
file of the singular isothermal sphere that we have already 
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used in paperl (n oc r~ 2 ), and which has not been much al- 
tered by the merging-process of the two BHs (Fig. 2). With 
n being the number density in the distance r = 1 pc we 
can write 



n(r)=no(y) . 



(11) 



Integrating this density from the inner to the outer radius 
of the torus (rj n and r out respectively) gives the surface 
density: 



/Tout 

/ n(r)dr 
Jr in 



n r 



r ut J 



To determine the number density no, we have to cal- 
culate the number of obscuring stars in the torus, which is 
the fraction 0.01 of the bound stars iV bn . Integrating the 
number density in Eq. (11) over the volume of the torus 
gives us the number of obscuring stars N os : 



27T rK-Otre /Tout 

sin#d# / n(r)r 2 dr 

J0 trB Jr ia 
2/ 



= 4?rcos6' trs nor (r ou t - r in ) 



where 9t rs is the half-opening angle of the torus. Solving 
this equation for no and applying it to Eq. (12), the wind 
radius is obtained with the help of Eq. (10), 



4 cos tTS r in r out 



1.8 x 10~ J pc. 



(14) 



For the parameters we assigned in the last step the val- 
ues for the case of a mass-ratio q = 1: The inner radius is 
r- m = 1 pc (see paper I and Sect. 2), and 5 pc is used for the 
outer radius since the density drops quickly with increas- 
ing distance. The half-opening angle is 60° (see Fig. 5) and 
iV os = 0.01iV bn , with N hn = 3 x 10 8 being a little more 
than the minimum of stars required to allow the BHs to 
merge (2.3 x 10 8 ). The wind radius computed in Eq. (14) is 
in very good agreement with our crude estimate in Eq. (9) 
and thus confirms the assumptions we made to obtain it. 

Such a wind has to be optically thick (Nn ~ 
10 24 cm~ 2 ) for all line of sights through its cross-section. 
We assume, that the density of such a tail, if projected 
along its radial extension, yields a constant surface den- 
sity. Therefore the mass in the wind can be estimated to 
be of order of 



M w = ttB^Nu to p « 8.4 x 10~ 2 M . 



(15) 



The final picture of our model for the torus is illus- 
trated (not to scale) in Fig. 8. The sum of the winds of 
the obscuring stars, shaped into elongated tails by the 
central radiation pressure, form a patchy, optically and 
geometrically thick torus. For an observer whose line of 
sight (LOS) is aligned with such a wind, the nucleus is ob- 
scured and only radiation scattered into the LOS will be 
detected. Since the number density of the stars decreases 
for smaller inclination angles between the LOS and the 
symmetry axis of the torus, the probability to observe the 



nucleus unobscured, or only partially covered, through a 
gap in the torus, increases. Such gaps will be continuously 
opend and closed due to the motion of the stars. According 
to our prerequisites Eq. (10) is satisfied and consequently 
the complete torus is optically thick on average. 

With the number of obscuring stars and the mass con- 
tained in their winds we can estimate the order of magni- 
tute of the dust-mass which is confined to the stellar winds 
outside the evaporation radius. The relation between the 
dust and gas mass is given by 



( 12 ) m d n d = lAZ d m p nn , 



(16) 



where m d and m p denote the mass of a dust grain and 
the proton mass respectively. The number-densities are 
accordingly given by n d and «h with Z d as the dust-to-gas 
mass-ratio. Integrating over the volume of the wind yields 
M d = 1.4Z d M gas . This is substituted in M w = M d + M gas 
and solved for the dust mass: 



(13) M d 



lAZ d 



1 + 1.4Z d 



-M, 



1.2 x nr 4 M ( 



© ■ 



(17) 



In the last expression we used the wind mass of Eq. (15) 
and for the dust-to-gas mass-ratio Z d w 10~ 3 (Hoefner 
et al. 1996). Multiplying this dust mass with the number 
of obscuring stars needed to maintain a covering factor of 
1 {N OB = 3x 10 6 ) gives the total dust mass contained in 
the wind of the obscuring stars, 



M d 

.total — 

M d N os = 350 M e . 



(18) 



This is about a factor of 10 more than the minimum val- 
ues given by Sanders et al. (1989), which are required to 
fit the spectra. Within spheres of the radii 1 pc and 10 pc 
they obtain for the minumum of dust amounts 0.2 Mq 
and 20 Mq respectively, using for the grains the parame- 
ters given in Mathis et al. (1977) and Biermann & Harwit 
(1980). If all grains within the given radius are larger than 
~ 1 «m, or are not directly exposed to UV radiation, the 
required dust mass is ~ 10 times higher, in good agree- 
ment with Eq. (18). Such shielding of the outer regions of 
the elongated winds by the inner parts is a natural conse- 
quence of our model and diminishes the UV irradiation of 
the outer parts of the wind. Also, due to a covering factor 
of about 1 it is likely that the wind of a star is covered 
by that of another one at smaller distances to the central 
source. 

The picture which emerges from our simulations and 
the results is the following: If a BBH is surrounded by 
a stellar population of comparable mass, the stars which 
are ejected in violent interactions with the BHs are able 
to carry away enough angular momentum of the binary 
so that it hardens till finally gravitational radiation dom- 
inates the shrinking process of the BHs (paperl). These 
eventually coalesce after about 10 7 yr, having spent most 
of the time in the range when the ejection of stars domi- 
nates the hardening. On the other hand for the same num- 
ber of initial stars there are enough left which are bound to 
the binary and form a torus-like distribution which peaks 
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Red Giants & their winds 
constituting a torus 





Fig. 8. This picture shows (not to scale) a sketch of the obscuring stars and their winds comprising a torus with the BBH in its 
center. In this cross-section through the torus, whose density peaks in about 3 pc distance, the obscuring stars are indicated by 
the small dots. They are surrounded by their winds, which are elongated by the central radiation into comet-like tails (shaded) 
with a lateral extent of about 0.004 pc. The central radiation source with the BBH is marked by the two BHs in the center. An 
observer with the LOS close to the edge of the torus might see the center through a gap directly, while the view of an observer 
with his LOS close to the equatorial plane of the BHs is probably blocked by the winds of at least one obscuring star. 



in about 3 pc distance. In this section we could show that 
the amount of obscuring stars in the remaining cluster is 
sufficiently large so that their winds achieve a total cover- 
ing factor of about 1. These winds are are optically thick. 
Thus the patchy torus as a total is geometrically as well 
as optically thick so that our model can indeed account 
for a dusty torus in the center of AGN, as is demanded by 
the unification model. 

4.1. Relativistic boosting 

Within the solid angle of a relativistc jet, the radiation is 
boosted and so the flux density of the approaching radia- 
tion, seen under an angle 9 to the velocity of the moving 
source, scales as (Rybicki & Lightman 1979; Longair 1981) 

S = D 3 ~ a S', 

with a being the spectral index and D the relativistic 
Doppler factor 



The primed quantities refer to the comoving frame of the 
source. Thus, for an index a ~ — 1, the observed flux den- 
sity depends on the Doppler factor to the power of four. 
Half of the radiation is emitted in a cone of half-opening 
angle 9 ~ I/7, if 7 ^> 1, and thus is strongly beamed. 
In this limit (7 > 1) we have D ss 7 and the luminosity 
in the cone is a strongly increasing function of 7. If we 
assume the luminosity L to scale in the same way as S 
with 7, it increases very strongly in the beaming cone (i.e. 
by a factor of 10 4 for 7 = 10 and a = —1). Consequently, 
according to Eq. (9), the extension of the winds of the ob- 
scuring stars, exposed to the strong beamed radiation, is 
negligible and the covering factor in this region tends to 
zero. This means that within the beaming cone, and if the 
jet is precessing, within the cone of precession, the emis- 
sion from the jet frees the polar cap regions from obscur- 
ing clouds, and therefore supports the toroidal structure 
of the central absorber. 



D 



1 



7(1-/3 cos 6>) 
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4.2. Lifetime of the torus 

After the merger of the BBH is completed, the toroidal 
structure of the stellar distribution will not collapse and 
is stable (paper I) . Therefore the torus will decay on 
timescales of the lifetime of its constituents, the obscuring 
stars (i.e. red super-giants and bloated super-giant stars) 
or has to be replenished with stars from outside. Because 
after the merger no torque acts anymore on the stars, the 
stars fed to the torus from larger distances would have 
to be accreted in an axi-symmetric way to maintain the 
torus-like shape of the circum-nuclear stellar distribution. 
The jet-outflow probably helps to keep the polar-cap re- 
gions free from stars. But probably the rate of stars ac- 
creted to the torus will decrease with time, as the outer 
regions, stirred up by the merger, will relax with time. 

Hence, unless the required number of obscuring stars 
is not maintained for a longer time by the evolution of the 
other stars in the cluster, the life-time of the torus will 
be of comparable order as the time needed by the BHs 
to coalesce from the distance where a torus emerged due 
to their torque acting on the surrounding stars. Once the 
binary has become hard (a <~ lpc), the BHs merge on 
scales of about 10 7 yr (paper I). 

In the next sections we present implications and pre- 
dictions of the proposed model. 

5. Jet-flip due to spin-flip of the primary BH 

As we have already mentioned in paper I, we expect the 
merger of two galaxies of comparable size to result in a 
galaxy of early type, i.e. of elliptical shape (e.g. Tinsley 
& Larson (1977), Biermann & Shapiro (1979), Farouki & 
Shapiro (1982)). The spiral structure the parent galaxies 
might have had will not survive such an event. On the 
other hand, radio-loud AGN so far have not been found 
in spiral host-galaxies, what suggests radio-loud galaxies 
to be the product of one or more major mergers. Since ob- 
servations show that spectra from AGN in the range from 
IR to X-rays, dominated by emission from an accretion 
disk, are basically the same (Sanders et al. 1989), Wilson 
& Colbert (1995) argue that the distinction between radio- 
loud and radio-quiet AGN is due to the different spins of 
the central supermassive BHs. They propose radio-loud 
objects to harbour a fast spinning super-massive BH as 
result of a recent major merger, whith the assumption of 
non- rotating BHs in the parent galaxies. The time scales 
we obtained for the merging of the BHs are considerably 
smaller than a Hubble-time, in agreement with such a 
model. 

But what happens, if one of the progenitor galaxies 
has a fast spinning BH in its center with a powerful radio- 
jet aligned with it's spin? The situation before such BHs 
merge is depicted in the left panel of Fig. 9. The spins of 
both BHs and the orbital angular momentum of the BBH 
are randomly oriented to each other, because the merger 
of both galaxies does not proceed in a favored plane rela- 
tive to that of the galaxies. While the secondary BHs spin 



is assumed to be negligible, we adopt for the primary BH 
a spin that is sufficiently large to power a strong radio-jet 
according to the spin-paradigm [L\ 3> £ 2 )- Its upper limit 
is given by MicR g = GM\jc, with R g being the gravita- 
tional radius. For the orbital angular momentum we have 
L m = Mxy/GM^/ y/q(l + q), where q = Mi/M 2 > 1 is 
the mass-ratio of both BHs and a is the semi-major axis. 
Thus the ratio of these two angular momenta is 



Li ^q{l + q) V R a 

It is a very interesting question now, what actually dom- 
inates the final spin of the merged black hole, L # . When 
the distance of both BHs has shrunk to about the last 
stable orbit (R\ s = 6i? g in case of zero spin, and in case 
of a Kerr BH it depends on the inclination, see Chirvasa 
(2002)) due to the emission of gravitational radiation, gen- 
eral relativistic effects become important. At this distance 
the orbital angular momentum dominates over the max- 
imum possible spin of Mi, up to a mass-ratio of q = 2. 
This is actually a conservative limit since we neglected 
the last stable orbit around M 2 . Hence the spin of the 
merged BH might be dominated by L m rather than L\, 
with the merged black hole M, spinning close to its max- 
imum value L, = GM, /c. This situation is shown in the 
right panel of Fig. 9. A new powerful radio-jet is emanat- 
ing from the center, aligned with the BH's spin which is 
pointing in the direction of L m , and therefore parallel to 
the symmetry-axis of the torus. This means that the jet 
flips into a new direction. The angle enclosed by the old 
and new jet-axis corresponds to that between L\ and L„ 
and hence can assume any value. After the merger the old 
jet is not fed anymore and its lobes will slowly fade away, 
while the new jet is digging its way through the ambient 
medium and the polar caps of the torus. 

This might be the explanation for some of the most pe- 
culiar objects seen on the sky, the X-shaped radio-galaxies. 
They show a secondary, presently non-active pair of lobes, 
which is larger than and sometimes almost perpendicular 
to the presently active pair (Parma et al. 1985; Rottmann 
et al. 1998; Rottmann 2001), which we interpret in this 
picture as the new jet. This scenario is also in agreement 
with the involved time scales of some 10 7 yr for the merger 
and - 6 x 10 7 yr in B2 0828+32 for the spectral aging of 
the secondary lobes ((Rottmann 2001), priv. comm. H. 
Rottmann) . 

Such a major merger will also increase the accretion- 
rate onto the BH which might be in favour of the argu- 
ments of Meier (2001), who suggests an association of jet 
production with geometrically thick accretion flows and 
BH rotation. 

For larger mass-ratios (q > 2) the possible influence 
of L m on the final BH spin is decreasing, and therefore 
also the bending of the jet into a new direction will have 
a smaller amplitude. Thus we do not expect a jump of 
the jet in minor mergers where M 2 <C Mi and therefore 
L„ <C Li (see Eq. (19)) and ascribe the phenomenon of 
X-shaped radio-galaxies to recent major mergers, where 
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Before merging: After merging: 




Fig. 9. This figure illustrates the change of the direction of the spin of the BH, induced by the merger of 2 massive BHs, and 
consequently the change of the direction of the jet. The left panel shows the situation before the merger, when the jet is aligned 
with the individual spin of the primary black hole of the binary system. The orbital angular momentum Z/„ of the BBH and the 
spins of both BHs (Li and L2) are randomly orientated to each other, since no direction of the plane of the merging galaxies is 
preferred relative to the spin of both the BHs in their centers. After the BHs coalesced (right panel) , the different spins combine 
to that of the merged BH with a new direction (L.). The jet emanating from the merged BH will be aligned with the new spin, 
which will be dominated by L s , and therefore has to jump form its orientation before the merger (Li) into this new direction. 



one of the progenitors already had a powerful radio-jet. 
If radio strength is correlated with the spin of the BH, 
then the observations imply that a BH merger of spinning 
black holes leaves the ratio of spin to mass J/M large. 

Before the two BHs finally coalesce, the surrounding 
patchy torus will have emerged on time scales of 10 yr 
(paper I) in the plane of the merger, with its symmetry 
axis pointing along the binary's angular momentum. As 
long as the semimajor axis of the BBH is much larger 
than the last stable orbit of Mi, the old jet of the pri- 
mary AGN is unaffected by M2 and still fed, if the feed- 
ing mechanism from the accretion disk is not interrupted. 
Because the orientation of this jet is not correlated with 
that of the torus, it might flow into the solid angle cov- 
ered by the torus, where it will interact with stellar winds 
and ISM in between. After the BHs coalesced the young 
past-merger jet is digging its way through the polar cap 
region of the torus, flowing along its symmetry axis. So 
both jets, the old and new, are interacting on sub-kpc 
scales with the ambient clumpy medium producing strong 
radio emission. This might be the cause of the so called 
Compact Symmetric Objects (CSO), which is a class of 
powerful radio sources consisting of high unbeamed lu- 
minosity radio emission regions separated by less than 
1 kpc and situated symmetrically about the center of ac- 
tivity (Phillips & Mutel 1982; Wilkinson et al. 1994). It 
is thought that the high-brightness regions are due to hot 
spots and mini- lobes, which are created by the termination 
of jets streaming out into opposite directions from the cen- 



ter, see Owsianik & Conway (1998) and references therein. 
Instantaneous speed variations of some components in 
CSO 0710+439 are most simply explained by interactions 
with a dense cloud medium, while other components are 
rapidly advancing through an intercloud medium. The ex- 
pansion velocities are all of order 0.2 h~ 1 c and therefore 
the age of these sources is estimated to be a few thou- 
sand years, showing that they are young rapidly grow- 
ing sources (Owsianik et al. 1998, 1999). It is concluded 
that CSOs are probably young extragalactic radio-sources 
which will evolve via Medium-size Symmetric Objects into 
Large-size Symmetric Objects such as lower luminosity 
FR II double radio-sources (Fanti et al. 1995; Readhead 
et al. 1996; Owsianik & Conway 1998). 

If the X-shaped radio galaxies are correctly explained 
with our proposed picture, the resulting gravitational 
wave pattern will strongly depend on the orientation of 
the three relevant spins (Chirvasa 2002) and therefore on 
the re-orientation of the primary black hole. It experiences 
a spin-flip in this picture. 

6. A patchy torus - explanation for BAL QSOs 
and transition from Type 1 to 2 

The broad absorption line (BAL) quasars comprise about 
10% of the optically selected quasars (i.e. Antonucci 
(2001b) or Schmidt & Hines (1999)), but since they are 
hard to find at optical wavelengths, their real fraction is 
thought to be in the range 20% — 30%. Moreover the BAL 
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covering factor is thought to be much less than unity and 
consequently there must be many objects, if not all, being 
intrinsically the same as BAL quasars, but not classified 
as such if seen from other directions. The region which 
is responsible for the formation of the BAL partially ab- 
sorbs the broad emission lines, and therefore it must be 
outside the BLR. Also the polarization in the broad emis- 
sion lines is observed to be lower than the polarization in 
the absorption troughs of the BALs. 

Inside the torus, with its inner edge close to the evap- 
oration radius of dust (<~ 1 pc) the heated matter will be 
accelerated and streams away from the nucleus through 
the gap between jet and torus cither in form of clouds or 
as a wind. Here, above the surface where the density of 
the torus gradually decreases, the clouds are exposed to 
more intense radiation. These arguments fit well into the 
unification model, where the BAL QSOs are quasars seen 
at intermediate inclination angles, with the line of sight 
grazing the surface of the obscuring torus (Weymann et al. 
1991; Voit et al. 1993) and thus is supported by our model. 

With the assumption that the fraction of BAL QSOs of 
X = 0.2 to 0.3 represents the covering factor Cbal in the 
same range (Weymann 1997), we can estimate the solid 
angle they comprise. We just have to solve the equation 
of the covering factors, \ = Cbal/CqsO; for the polar 
angle 6>bal , using the relations Cqso = 1 — cos #trs and 
Cbal = cos 6*bal — cos 6> trs • In this notation the QSO ex- 
hibits BAL features if seen under inclination angles in the 
range #bal < #inci ^ #trs- Assuming the half-opening an- 
gle of the torus to be 6>t rs = 45° and using the above frac- 
tions of BAL QSOs x = 0.2 and 0.3 we obtain 6> B al = 40° 
and 37° respectively. Hence for an intermediate inclination 
angle in the range of about 

38° < #inci < 45° 

the observed quasar will appear as a BAL QSO, in agree- 
ment with Schmidt & Hines (1999) who find BAL QSOs 
to be typically 2.4 times more polarized in the optical than 
QSOs. 

Among the BAL QSOs the subclass of low-ionization 
absorbers shows evidence to be stronger polarized and to 
be reddened by dust (Sprayberry & Foltz 1992; Egami 
et al. 1996). Green et al. (2001) assumed that the same 
intrinsic power-law with index F w 1.8, which is con- 
sistent with the mean slope of radio quiet QSOs and 
has been derived for high-ionization BAL QSOs, ap- 
plies to all BAL QSOs and is partially covered. This 
assumption has been confirmed in some recent publica- 
tions (Gallagher et al. 2001; Green et al. 2001; Gallagher 
et al. 2002). For low-ionization BAL QSOs, which are 
weak in X-rays, Green et al. (2001) infer that they are 
enshrouded by an additional intrinsic column density of 
nearly 10 23 ~ 24 cm~ 2 . These low-ionization BAL QSOs 
have been proposed by Brotherton et al. (1997) as the 
most edge-on QSOs. 

For larger inclinations than the intermediate range the 
line of sight will intersect with more absorbing material 
since it approaches the surface of the patchy torus, which 



does not have a sharp defined edge. Consequently less of 
the BLR is seen directly and the fraction of scattered light 
as well as the polarization is increasing. With increasing 
inclination the view to the nucleus is increasingly blocked 
by the torus and the quasar will appear more reddened 
(IRAS QSO) till finally, if seen edge on, it will appear as 
a QSO of Type 2 (ULIRG or FR II radio galaxy). On the 
other hand, objects seen at lower inclinations don't suffer 
absorption of the central radiation by the torus and are 
classified as quasars. 

In the special case of the quasar FIRST 
J101614.3+520916 Gregg et al. (2000) find BAL 
features in the optical and radio loud emission from lobes 
of classic FR II type at the same time. According to 
the authors this quasar has some properties in common 
with low-ionization BAL QSOs, as for example higher 
reddening. They deduce an inclination of the jet-axis 
to the LOS of more than 40° and conclude that this 
quasar is in contradiction with the orientation model as 
explanation for BAL QSOs. Gregg et al. (2000) suggest 
that this quasar rather happens to be seen in a rare 
and short-lived state, showing both, BAL features and 
developed radio lobes at the same time, and postula- 
tion that it is a rejuvenated quasar, possibly through 
merger. Contrary to this conclusion we claim that these 
observations actually support the orientation model. 
According to the above numbers, BAL QSOs are seen 
at intermediate inclination angles in the range between 
~ 38° and <~ 45°, in very good agreement with the 
inclination obtained for this quasar. Thus the LOS to the 
core of J101614. 3+520916 is grazing the surface of the 
patchy torus, whose constituting absorbers, the stellar 
winds, arc increasing in number density towards the 
equatorial plane. At higher inclinations the torus would 
block the free view to the center completely and no BALs 
could be detected, so that this object would appear as 
a typical FR II quasar. This is also in agreement with 
BAL QSOs, that show large radio lobes, being so rare. 
We do not exclude evolutionary effects, but as shown 
above, we also expect a merger to seriously affect the 
jet, what should be detectable. And the question to be 
answered for the evolutionary/rejuvenated scenario is: 
What is the explanation then for the difference between 
the radio-quiet BAL and non-BAL quasars? 

Due to the patchy nature of the torus, we propose par- 
tial covering absorption plus scattering in order to fit the 
spectra of BAL QSOs. Another important consequence of 
the patchy torus is, when the absorbing clouds are moving 
through the line of sight, the column density changes and 
variations in the spectrum are expected. This is then a 
natural explanation for the variability of the BAL quasar 
PG 2112+059 between the ROSAT and ASCA observa- 
tions on scales of 6yr, being due to changes in the ab- 
sorber, either the ionization parameter, or the column den- 
sity. We investigate such variabilities in more detail in the 
following. 

One observable consequence of an obscuring patchy 
torus, as developed in this paper, is a possible variation 
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of the absorption properties, i.e. the column density. A 
covering factor of the torus of about unity means that 
not necessarily all lines of sight through the torus to the 
central radiation source are covered and that due to the 
motions of the stars in the torus under the influence of 
the central BHs such gaps will form and will be closed 
continuously. 

When an obscuring cloud or stellar wind of the torus 
is moving into (out of) the line of sight of the observer to 
the radiation source, the column density will drastically 
increase (decrease) and the source will gradually appear 
much weaker (brighter) in luminosity. For our chosen AGN 
with Mi = 10 8 M© we obtained the circumnuclear torus 
to be located in a distance of a few parsecs to the center. 
The radius of the wind of the obscuring stars we computed 
to be of order of 0.002 pc. For a star moving on Keplerian 
or bits in t he potential of such a BH, the velocity is v = 
yjGMi/r. Therefore the time required for the star's wind 
to move along its diameter from one edge to the other 
through the line of sight is 



. Rv. 



10 



(JL\ 



0.002 pcV3pcy V 1Q8M e 



Mi 



yr. 



(20) 



Consequently variations of the column density and the lu- 
minosity are expected to happen on scales of a decade. 
Since the density of stars is enhanced towards the equa- 
torial plane of the torus it is likely that there is more 
than one star in the line of sight and thus a transition 
from optically thick to thin or vice versa will be a rare 
event in an edge-on seen torus. But if the LOS grazes the 
surface the torus the chance to observe such a transition 
will be higher. And indeed there are a couple of observa- 
tions which fit well to this interpretation. Before compar- 
ing the above computed time scale for a transformation 
from Type 1 — » 2 with observations, we want to scale it 
just in a simple way with the mass of the central BH. 

Assuming the luminosity to be proportional to the ed- 
dington, L c dd = 47rGcMim p /<7Th, we have L oc Mi. Thus 
with Eq. (9) the wind radius scales as 



Together with v k op oc y/Mi/r the dependency of the time 
on mass and radius is: 



oc 



r 5/2 
3/2 ' 



Mi 

There are now at least two possibilities for the choice of 
the inner radius of the torus: 

(1) The evaporation radius, for which we have r evap ex 
L 1 ! 2 oc M\ 12 (sec Eq. (7)). Thus we obtain the depen- 
dency 



what means that the time required for a cloud passing 
through the line of sight is decreasing with increasing mass 
of the central BH. 

(2) Another choice for the radius of the torus is the 
distance of the BHs, when the binary is hard, since the 
simulation showed that this coincides to within a factor of 
a few with the radius of the torus (paper I) . According to 
Milosavljevic & Merritt (2001) the binary becomes hard 
at a semi-major axis ah = GM m j%a 2 . Using the relation 



M. = 1.3 x 10 8 M G 



200km/s 



with a = 4.72(±0.36) (found by Meritt & Ferrarese 2001a; 
2001b) we get 
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Again there is a negative exponent, albeit small, so that t 
is decreasing with increasing Mi , as for r — r ovap . 

Since the exponents in the time-mass relation are much 
smaller than 1 for both choices of the inner radius of the 
torus, the transition-time does not depend strongly on the 
mass of the central BH, and is almost constant on orders 
of 10 yr. For a patchy absorber it is possible that more 
than one cloud or stellar wind happen to block the line 
of sight, and depending on their individual velocities and 
directions in which they move, it is possible that the de- 
crease and increase of the central flux proceeds on different 
time scales. 

In the following we will give some examples of sources 
which are in agreement with our patchy torus model 
within the unification scheme and show the expected vari- 
ations in column density on scales of a decade. 

NGC 7582: Compared to the EXOSAT observations 
in 1983 Warwick et al. (1993) detect 4yr later with 
Ginga a significant increase of the column density to 
~ 4.6 x 10 23 cm~ 2 by a factor of about 3. While the 
spectral index between the two ASCA observations from 
November 1994 and 1996 does not show a significant varia- 
tion the column density seems to have increased by ~ 44% 
(Xue et al. 1998). In the hard and soft X-ray band of the 
BeppoSAX observation in November 1998 Turner et al. 
(2000) find the nucleus clearly to have brightened since 
the 1994 ASCA observation. While this increase seems 
to be consistent with the appearance of holes in the full 
screen (1.4 x 10 23 cm -2 ) the authors also find evidence for 
larger absorption at the BeppoSAX epoch with a "thick 
absorber" covering 60% of the nucleus, corresponding to 
a half-opening angle of a torus of ~ 53° , which has a col- 
umn density of 1.6 x 10 24 cm~ 2 . Both these values are in 
good agreement with the column density we computed for 
the torus and that what is obtained if the matter would 
be smoothly distributed in the volume of the torus (see 
Eq. (23) in Sect. 7). The total timescale of about 15 yr 
for major variations in the column density is in very good 
agreement with our predictions. 

NGC 2992: This source has been monitored over a 
time of 20 yr with different telescopes (HE AO (Mushotzky 
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1982; Singh et al. 1985), Einstein (Halpcrn 1982; Turner 
et al. 1991) EXOSAT, (Turner & Pounds 1989), Gmga 
(Nandra & Pounds 1994), ASCA (Weaver et al. 1996) and 
BeppoSAX (Gilli et al. 2000)). During 16yr the 2-10 keV 
flux declines steadily by a factor of 20 and is 4yr later 
back to its initial value measured in 1978. Over this time 
range the spectrum retains a constant shape, which is only 
different in the low state observed by ASCA in 1994, when 
there arc also no short-time variabilities are seen. This is 
actually what is expected, since according to our idea in 
this stage the direct view to the center is blocked by a 
cloud. 

NGC 4051: While the flux in the range 2 - lOkeV of 
this Seyfert 1 galaxy has been almost constant for more 
than 10 yr till the RXTE observations by Uttley et al. 
(1998) it is measured to be 20 times fainter 1.5 yr later 
with BeppoSAX in 1998 by Guainazzi et al. (1998). The 
time scale of ~ 1.5 yr for the decrease of the luminosity by 
a factor of 20 might be too short to be accounted for by 
a cloud of the patchy torus moving into the line of sight. 

NGC 1365: Between the ASCA observations and the 
BeppoSAX in August 1997 three years later the flux in the 
range 2—10 keV has increased by a factor of 6 (Risaliti 
et al. 2000). 

NGC 3227: In their comparision of the data from 
ASCA observations in 1993 and 1995 and ROSAT obser- 
vations in 1993 George et al. (1998) show that the column 
density of the ionized absorber increased by about an or- 
der of magnitude, i.e. from Nh ~ 3x 10 21 cm -2 in 1993 to 
]\f H ~ 3 x 10 22 cm" 2 in 1995. They think this to be most 
naturally explained by a cloud of material moving into the 
cylinder of sight. 

7. Influence and implications of the binary's 
mass-ratio 

Our simulations in paper I showed that for a decreas- 
ing mass-ratio of the BHs (q = M\jM<i > 1, with 
Mi = 10 8 Mq = const.) the fraction of the ejected stars 
is increasing. But to extract the bigger amount of an- 
gular momentum from a more massive binary (L^ = 
a/ GMfa/ \J q(q + 1)), the increased fraction of ejected 
stars is not sufficient, so that the star cluster has to be 
more massive in order to absorb the binary's angular mo- 
mentum and to enable the BHs to coalesce. Therefore the 
amount of bound stars has to increase with the mass M 2 
of the secondary BH, and so the torus becomes more mas- 
sive. 

As we have pointed out in the discussion of the influ- 
ence of the mass-ratio, the inner regions become increas- 
ingly unstable for more massive secondary BHs, so that 
at the inner edge of the torus the stars are moving on al- 
most circular orbits to avoid violent interactions with the 
binary. For sufficiently small mass-ratios a more sharp de- 
fined torus forms during the late stages of the merger, 
while for young mergers and mergers with large mass- 
ratios a more diffuse and shell-like density distribution of 
the stars is maintained. This is clearly visible in Figs. 2 



and 14 of paper I, which also show a larger half-opening 
angle of the torus in the range 50° to 60° for q = 1 (i.e. 
M 2 = Mi) compared to about 45° for q = 10. This is due 
to the stronger torque that is exerted by the binary with 
smaller q on stars in the polar cap region, and therefore 
less orbits are passing through this region. 

In major mergers we also expect the accretion rate on 
the central BH to be much higher and thus these merg- 
ers to be much more luminous than minor ones. This is 
in agreement with the idea that the opening angle de- 
pends on the luminosity, being larger for the more lumi- 
nous AGN. Consequently the radiation pressure acting on 
the torus and its stars is less in minor mergers, and accord- 
ing to Eq. (9) the radius of the stellar winds can extend to 
much larger radii, not being stretched into elongated tails. 
Hence the column density along the LOS to the center is 
diminished and the torus becomes optically thin (see also 
Sect. 4). 

This can be illustrated with a simple order-of- 
magnitude estimate. The volume of the torus with half- 
opening angle 9 tIS within the radial limits r- m and r out 
is Vt rs = ^-{r^ nt — rf D ) cos 9 tIS - If we distribute the total 
mass contained in the stellar winds homogeneously in this 
volume, the number density of Hydrogene is 

M^total 

lAZ d m p V trs ' 

Integrating along the line of sight from r; n to r out , yields 
for the column density of such a homogeneous torus 

Nk = 3 M d , to tai r^t - r in = 4 g x 1q23 cm - 2 
5.6 7rZ d TOp cos fltrs r* ut - rf n 

(23) 

The total dust mass of the torus is 350 M Q (see Eq. (18)) 
and for the geometry of the torus we used the same val- 
ues as in Sect. 4, i.e. 60° for the half-opening angle and 1 
and 5 pc for the inner and outer radius respectively. This 
gives less than half the column density of the patchy torus 
in Sect. 4, when the radiation pressure is strong enough 
to turn back the wind into elongated tails. Thus a torus 
with smoothly distributed matter would be optically thin, 
showing that a sufficiently high luminosity is necessary in 
order to maintain an opaque torus. According to our re- 
sults from paper I the optical thickness of tori as products 
from major mergers (Mi = M2 = 10 8 M Q ) compared to 
minor mergers (Mi = 10 8 M Q , M 2 = 10 7 or 10 6 M Q ) is 
further enhanced due to the higher absolute amount of 
stars that stay bound in the potential of the binary, even 
though a larger fraction is ejected. This is plausible, since 
in a minor merger a less massive stellar cluster is surround- 
ing the secondary BH, which is dragged to the common 
center by dynamical friction. 

Thus we come to the following conclusion: In a ma- 
jor merger with two BHs of comparable mass (<~ 10 8 M Q ) 
more stars are brought to the common center and stay 
bound in the potential of the binary, where they en- 
hance the optical thickness of the torus they constitute. 
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Because of the stronger torque exerted by the BHs, the 
half-opening angle of this torus is larger than in a minor 
merger. Since the accretion rate is expected to be close 
to the Eddington limit, the stellar winds are exposed to 
strong radiation pressure from the center, which shapes 
the winds into elongated cometary tails, pointing radially 
away from the center. Hence they are sufficiently opaque 
as to make the stellar torus optically thick. With increas- 
ing time the accretion rate and consequently the luminos- 
ity will decrease, so that the stellar winds in the torus will 
become less collimated. Therefore the opacity of the torus 
will decrease also as an evolutionary effect of the AGN. 

If, on the other hand, a 10 s Mq BH merges with a black 
hole of much lower mass, fewer stars are involved and the 
torus will contain less stars. Because of the smaller torque 
of the binary, the torus is more diffuse and the opening 
angle is smaller. According to our simulation, such a torus 
looks like that of a major merger in early stages (compare 
the Figs. 2 and 14 in paper I). In a minor merger also 
the accretion rate will be smaller and consequently the 
luminosity. This results in less radiation pressure acting 
on the stellar winds, which are not shaped into elongated 
tails, and thus the opacity of the torus, containing already 
fewer stars, is even more diminished. 

This seems to be in agreement with the observations: 
Major mergers of galaxies of comparable size with BHs in 
their centers of comparable mass will violently stir up the 
stars, gas and dust, and finally assume an elliptical rather 
than a spiral shape. According to our finding above, ellipti- 
cals then should harbour AGN with tori having larger half- 
opening angles than spirals and thus on average Type 1 
nuclei should be detected more often in elliptical than spi- 
ral hosts. This is supported by the survey of Malkan et al. 
(1998), who find Seyfcrt 1 nuclei to reside on average in 
more early type hosts than Seyfert 2 nuclei. 

It has been emphasized (i.e. Antonucci (2001a) and 
references therein) that there must be a range of covering 
factors for dusty tori, with the Type 2-classified objects 
having higher average covering factors. It is concluded that 
the populations are therefore intrinsically different in their 
statistical properties to some extend. Maybe the mass- 
ratio of the BHs and the evolution of AGN are the reasons 
for statistically different covering factors. 

If the nucleus happens to be oriented to us in a way 
that the line of sight grazes the edge of the torus, the 
probability to see the nucleus still directly through a gap 
in the edge of the clumpy structure is more probable for 
large than small mass-ratios. This might be an explanation 
for NGC 4151, which is a Type 1 and has an aligned NLR 
too. But HST images show a cone on subarcsec scales 
(Antonucci 1993). According to the strict unified model 
this is not expected since a Type 1 should only be seen if 
the observer is inside the unobscured solid angle so that 
a projected ionization cone can not be seen. In our model 
the patchy torus can have gaps or holes, the more likely 
the larger q is and the closer the line of sight comes to 
the edge of the torus which is not clearly defined. In the 
case of NGC 4151 the observer seems to be outside the 



opening angle. But through a gap in the torus close to its 
edge enough radiation can escape and is seen directly so 
that NGC 4151 is classified as a Type 1. 

M87, the nearest giant elliptical galaxy (z = 0.0043), 
has been observed at 10.8 /im wavelength by Perlman et al. 
(2001). If there is a dusty opaque torus it should be vis- 
ible in the reradiated IR, but the authors find only little 
evidence of thermal emission from dust. Because of its 
large BH mass (~3x 10 9 M Q , Marconi et al. (1997)) and 
the giant elliptical shape it is likely to have undergone a 
major merger in the past. Since then the fuelling of the 
AGN probably has decreased and the luminosity dropped 
to its current estimated value of 10 42 erg/s (Whysong & 
Antonucci 2001), what fits to M87 being on the FR I- 
II border, as well in morphology as in radio power (Owen 
et al. 2000). This luminosity is too weak as to maintain an 
opaque torus, which also might have been dissolved since 
the merger (see Sect. 4.2, lifetime of the torus). A possi- 
ble minor merger afterwards most likely would not have 
sufficiently increased the luminosity and replenished the 
torus with stars in order to form a dusty and opaque torus 
again. Whysong & Antonucci (2001) measured the rera- 
diation from nuclear hot and warm dust, which must be 
emitted according to the unification scheme by the torus 
that absorbs the radiation from the AGN it surrounds. 
This radiation should be found in almost any hidden AGN 
and provides an estimate of the unblocked central lumi- 
nosity. While for Cyg A at rcdshift 0.057 (more than 10 
times the redshift of M87), with an estimated luminosity 
~ 1.5 x 10 45 erg/s, they could detect this component, it 
is much weaker in M87, where no hidden nucleus can be 
found. The estimated mass accretion in Eddington units 
M/M c dd is in the range 10~ 35 to 10~ 3 , which is much 
less than 1/50 (Liu et al. 1999). Below this value a thin 
Shakura-Sunyaev disk can not survive (Meyer et al. 2000). 
With respect to its large mass, M87 seems to host a starv- 
ing black hole in the center. 

This is in line with the results of Meisenheimer et al. 
(2001), who compared galaxy-quasar pairs from the 3CR 
catalogue in the IR range from 5 to 180 /im to test the uni- 
fication scheme for luminous radio galaxies and quasars. 
The pairs have been selected such that they match in 
178 MHz luminosity, which is thought to be emitted fairly 
isotropic, and rcdshift in order to minimize the effects 
of cosmic evolution. They can not distinguish the pairs 
by their mid- and far-infrared properties, what strongly 
supports the unification scheme. The authors also find 
the ratio of thermal dust power vF Vl averaged over 60 
and 100 fim, to the radio power at 178 MHz to correlate 
better with redshift than with luminosity. They suggest 
that this might be due to the thermal power of a radio 
source being primarily controlled by the accretion rate, 
with sources accreting at high rates being more numerous 
at large redshifts, when major mergers have been more 
frequently than today. 
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8. Summary and conclusions 

In paper I we showed that a central geometrically thick 
torus, comprised of stars, results at a distance of order 
3 pc from the center as a product from the merger of two 
galaxies and their central supermassive black holes. This 
stellar torus has to be about as massive as the binary black 
hole in order to enable the BHs to get rid of their orbital 
angular momentum and to coalesce on scales of 10 7 yr. 

In the present article we proved that this torus with its 
patchy structure is in very good agreement with the prop- 
erties of the ubiquituous torus in AGN, as are demanded 
in the unification scheme and deduced from observations. 
About 1% of these stars have strong enough winds in or- 
der to obscure the central radiation source for lines of 
sight passing through their winds. The central radiation 
pressure shapes the winds of these obscuring stars in the 
torus into elongated tails, pointing radially away from the 
central source. Consequently their azimuthal extension is 
about 4 x 10~ 3 pc, such that these winds along their tails 
are optically thick, and that their total covering factor 
within the torus amounts to ~ 1. Thus these winds cause 
the geometrically thick torus also to be optically thick 
with column densities of about 10 24 cm -2 , just as is ob- 
served. Also the inner radius of about 1 pc coincides with 
the evaporation radius of graphite dust grains (Lawrence 
1991), what has been previously assumed to be the inner 
edge, and what is also the distance where the BHs become 
hard (Milosavljevic & Merritt 2001) and which defined the 
inner edge in the simulation of paper I. 

In Sect. 5 we demonstrated that the recent merging 
of two comparable supermassive BHs, the prerequisite for 
the forming of the opaque torus, can explain the X-shaped 
radio galaxies. According to the spin-paradigm, fast spin- 
ning BHs are powering strong radio jets in radio galax- 
ies. If such a galaxy merges with another one hosting a 
BH of comparable mass, the orbital angular momentum 
of the resulting BBH dominates over the maximum possi- 
ble spin of the primary BH, at least till their major axis 
has shrunk to the last stable orbit. When the BHs even- 
tually merge, the final spin might be dominated by the 
orbital angular momentum, leaving the merged BH spin- 
ning close to its maximum value. As a consequence a new 
powerful jet emanates into the direction of the orbital an- 
gular momentum, which is aligned with the symmetry axis 
of the torus. Therefore this jet streams through the am- 
bient medium, producing strong radio emission, while the 
old jet's lobes, not fed anymore, are slowly fading away. 
The time scales of spectral aging of these lobes is close to 
the merger-time of the BHs and thus supports further this 
interpretation. The old jet might intersect with the torus 
and consequently would interact with its clumpy medium, 
also resulting in strong radio emission. This could be an 
explanation for the Compact Symmetric Objects, which 
are thought to be young and probably evolve into large- 
size Symmetric Objects, such as lower luminosity FR II 
double radio sources. 



With such a patchy torus-model, the BAL QSOs fit 
well into the sequence in which a QSO appears as Blazar 
if seen pole-on. As the inclination angle increases it ap- 
pears as a normal quasar, then high ionization BAL QSO, 
low ionization BAL QSO and finally, when the line of sight 
lies in the equatorial plane of the torus, as ULIRG. For 
the intermediate agles the line of sight grazes the surface 
of the torus, explaining the features of the BAL QSOs. In 
the sample studied by Gallagher et al. (2002), for most of 
them a partial-covering absorber provides a significantly 
better fit than other models, confirming our torus model. 
The detected variability of the BAL quasar PG 2112+059 
on scales of 6yr is ascribed to changes in the absorber, 
probably the column density. This is just the time the 
stellar winds in the torus need to move through the line of 
sight. In Sect. 6 we give more examples of objects whose 
column densities are observed to change considerably on 
these time scales. They are very strongly supporting the 
idea of the patchy torus, since the winds of the stars, 
moving in the potential of the central BH, have the right 
size and optical depth at the appropriate distance to the 
center to yield such strong variations in A^h on the right 
time scales, as they move through the line of sight. Due to 
the stellar motion gaps in the torus will be continuously 
opened and closed. 

For major mergers the resulting torus is more massive 
and therefore has a larger column density (see Sect. 7). 
Because the accretion-rate and consequently the luminos- 
ity will be higher than in minor mergers, the radiation 
pressure acting on the winds is stronger, being able to 
bend the winds into elongated tails. This further increases 
the torus' opacity compared to that of minor mergers. 
Since the torque of a more massive binary acting on the 
stars is stronger, the opening angle of the resulting torus is 
wider and therefore correlates with the central lumnosity, 
as has been suggested in the past. But the column density 
of the torus will also change with time, as the luminosity 
decreases when the accretion rate weakens. Thus the opac- 
ity depends on the mass-ratio of the merging black holes 
and on evolutionary effects of the AGN. Hence we ex- 
pect to observe Type 1 AGN more likely in elliptical host 
galaxies at larger rcdshifts, where major mergers occured 
more frequently We also expect higher covering factors of 
Type 2 objects. In their survey of Seyfert galaxies Malkan 
et al. (1998) find on average Type 1 AGN to be hosted 
by more early type galaxies than Seyfert 2s, in agreement 
with our reasoning. 

Previous compact torus models with a smooth distri- 
bution of the matter all faced the same problem in pre- 
dicting a too narrow infrared spectrum, and to fit the 
data, additional NIR sources had to be invoked (Pier 
& Krolik 1993; Granato & Danese 1994; Efstathiou & 
Rowan-Robinson 1995; Alonso-Herrero et al. 2001). As has 
been already mentioned in these papers, it is obvious that 
a patchy torus will tend to increase the dust temperature 
in the outer parts, since they are exposed to the radiation 
from the center shining through gaps in the inner parts of 
the torus. This will tend to broaden the IR spectrum. Such 



18 



C. Zier and P. L. Biermann: Binary black holes and tori in AGN 



models also exhibit strong emission or absorption 9.7 
Silicate features, which are usually not in agreement with 
the observations. Another problem is to achieve the geo- 
metrical thickness of the torus, which according to Pier 
& Krolik (1992b) can be supported by radiation pressure. 
But for a low luminosity source they need the torus to be 
clumpy and introduce a radiation pressure driven random 
motion of the clumps in order to maintain the thickness 
of the torus. 

In the model we proposed here and in paper I, the torus 
is naturally thick as a result of the stars moving in the po- 
tential of the binary. The patchy structure inevitably leads 
to self-shielding of the stellar winds and also within the 
winds. This probably has an effect on the temperature of 
the gas and dust as a function of the distance to the cen- 
tral source. Both, the selfshiclding and the temperature 
distribution in the torus will have a strong influence on 
the reprocessed and reemitted radiation. The calculation 
of the temperature distribution and spectra of the winds 
and the torus they constitute would involve a full three- 
dimensional treatment, so that radiation-transfer calcu- 
lations are well beyond the scope of this paper and we 
can not give a comment on possible Silicate features as 
function of the inclination angle. 

The torus model proposed in these two articles gives a 
coherent picture with respect to the formation of the torus, 
its evolution and application to observations in terms of 
the unification scheme. 
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